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We address the quantum statistics of electrons created in the low-energy edge-state Hilbert space 
sector of incompressible fractional quantum Hall states, considering the possibility that they may not 
satisfy Fermi statistics. We argue that this property is not a priori obvious, and present numerical 
evidence based on finite-size exact-diagonalization calculations that it does not hold in general. We 
discuss different possible forms for the expression for the electron creation operator in terms of edge 
boson fields and show that none are consistent with our numerical results on finite-size v = 2/^ states 
with short-range electron-electron interactions. Finally, we discuss the current body of experimental 
results on tunneling into quantum Hall edges in the context of this result. 



I. INTRODUCTION 



The quantum Hall (QH) eSec^^ occurs when an in- 
compressibility, i.e., a discontinuity in the dependence 
of chemical potential fx on density, occurs at a two- 
dimensional (2D) electron system (ES) sheet density tiq 
that is magnetic-field B dependent. The QH effect at in- 
teger Landau-level filling factors v — 2n£'^no arises from 
the quantization of 2D ES kinetic energy and from the 
macroscopic degeneracy of Landau-level states with a 
particular kinetic energy. {£ = \J'hcl\eB\ is known as 
the magnetic length.) The fractional QH effect (chem- 
ical potential jumps at fractional values of v)^ on the 
other hand, does not occur in a non-interacting electron 
system, and is due to constraints on the correlations^ 
that can be achieved among electrons that have the 
same quantized kinetic energy. A necessary consequenceSi 
of magnetic-field dependence in no is the existence of 
states in the chemical-potential gap that are localized 
at the edge of a finite-size system and carry equilib- 
rium current / with the property dJ/d/i = vejhA These 
edge-electron systems are obviously one-dimensional and, 
since they carry an equilibrium current, obviously chiraL^ 
Microscopically^, the edge of a non-interacting electron 
system at integer filling factor = to is equivalent, at 
low energies, to a one-dimensional electron system with 
TO flavors of fermions that can travel only in one di- 
rection, i.e., TO chiral fermion branches. It was argued 
some time a%o^^ on the basis of trial wavefunctions for 
finite-size systems, that the edges of incompressible frac- 
tional QH states can also be described microscopically 
as chiral one-dimensional electron fluids that have, in 
general, an unequal number of incquivalent left-moving 
and right-moving branches. In a series of beautiful pa- 
pers based on hydrodynamic and fleld-theoretic consid- 
erations, We n^i^i^" proposed and developed the idea that 
the properties of fractional QH edges could be described 



using a generalization of the bosonization approach that 
was developed earlier for conventional one-dimensional 
electron systemsiiti^iiii^ and provides a simple descrip- 
tion of their characteristioi^ Luttinger-liquid power-law 
correlation functions. In Wen's theory, edge excitations 
can be described microscopically in terms of a number 
of chiral boson fields, and the resulting edge system is a 
chiral Luttinger liquid (xLL) . Finite-size numerical calcu- 
lationsi^ have been a useful tool in verifying the funda- 
mentally bosonic character of the edge-excitation spec- 
trum and in testing some experimental predictions of 
xLL theory. Further experimental evidence in support of 
some aspects of xLL theory is summarized below. There 
are, however, difficulties in reconciling this effort to cap- 
ture generic aspects of the microscopic physics of frac- 
tional QH edges with experimental observations. Early 
work-^ did appear to imply that the edge structure of 
QH systems at the Laughlin series of filling factors, i.e., 
for V = l/(2p -|- 1) with positive integer p, is rather well- 
described in terms of a single-branch xLL characterized 
by a power-law exponent a = ^/v, as expected on the 
basis of xLL considerations. Recently, however, this uni- 
versal dependence of a on filling factor has been ques- 
tioned both theoreticalljiiSiiSiSfliSi and experimentallyi^S 
In addition, tunneling density-of-states observations^^ at 
hierarchical filling factors, i.e., v = ml{2mp± 1) with 
integer to > 1 are in apparent contradiction with the 
predictions24 of xLL theory. 

Motivated by this stark experimental discrepancy, we 
re-examine in this paper a key ansatz of X-^L theory, 
which has appeared to us to be non-obvious^ and con- 
cerns the properties of the operator t/i^ obtained from the 
full microscopic electron creation operator -0^ by project- 
ing onto the low-energy sector of edge excitations: 

^t = 7?^t7? ^ (1) 

where V = J2{rji} I {Vi}} {{Vi} I is the projection operator 
onto the Fock-space subset of low-energy (edge) excita- 
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tions I {iji}) , discussed at greater length below. In apply- 
ing xLL theory to evaluate electronic correlation func- 
tions, the representation of ip in terms of bosonic edge- 
density fluctuations is a key ingredient. Such bosoniza- 
tion identities can be derived constructivelji2SiSL2& for 
a conventional one-dimensional electron system^^ Their 
generalization to the fractional-QH casej2iS£ in which 
there is no adiabatic connection to non-interacting elec- 
tron states, must however be based on appealing but 
heuristic physical arguments that, ultimately, have to 
be verified by experiment, or by numerical calculations. 
In our view it is not clear beyond any doubt that the 
bosonized forms of ip that have been used in xLL the- 
ory to evaluate electronic correlation functions and pre- 
dict observables, like the power-law exponents in tun- 
neling current- voltage characteristics, are always correct. 
In particular, an important guiding principle that has 
been used to limit possible bosonized expressions for elec- 
tron field operators in xLL theory is the seemingly ob- 
vious requirement that they satisfy Fermi statistics. In 
this paper, we use finite-size exact diagonalization stud- 
ies of a short-range-interaction model to directly test the 
Fermi-statistics ansatz. The subset of Fock-space states 
that represent edge excitations of this model can, for the 
most part, be identified convincingly. We demonstrate by 
explicit calculation of some anticommutator matrix ele- 
ments that electron creation operators projected onto this 
low- energy Fock space do not satisfy Fermi anticommu- 
tation rules. In Section ^ of this paper, the numerical 
calculations that support this claim are described in de- 
tail. Section ITTll discusses the problem of understanding 
the properties of the projected electron creation operator 
and of finding a useful expression for it in terms of edge 
boson fields, in light of our numerical finding. We con- 
clude in Section llVl with a brief summary. A preliminary 
report on this work was presented earlieri^ 

II. NUMERICAL TEST OF THE FERMI 
STATISTICS ANSATZ 

The second-quantized operator [ilj{x,y)] that 

creates [annihilates] 2D electrons in the lowest Landau 
level obeys anticommutation relations that encode the 
fundamental antisymmetry condition satisfied by many- 
fermion wavefunctions: 

{^PHx,y),^p\x',y')}^0 . (2) 

The question we address in this section is whether the an- 
ticommutation relations are still satisfied after projection 
onto the low-energy (long-wavelength) sectors of Fock- 
space that represent edge excitations of particular incom- 
pressible states. Because of the projection, Eq. 101 does 
not mathematically guarantee the relation 

4,\x,y)4>\x',y') = -^/j\x',y')4'Hx:y) ■ (3) 

A physical argument along the following lines does, how- 
ever, appear plausible. It is possible to add two different 



electrons to the system at low energies that are localized 
at different positions along the edge via processes that 
can be represented by the projected creation operator. 
The edge-electron positions can then be adiabatically in- 
terchanged, ending up with an equivalent many-fermion 
state that must differ only by a sign from the original 
state. If the many-particle state can be represented, at 
all intermediate relative positions, by two projected cre- 
ation operators acting on the starting state, it seems hard 
to escape the conclusion that these operators must satisfy 
Fermi statistics. However, the correlations that establish 
the bulk gap could be disturbed when the two electrons 
are in close proximity. It is therefore difficult to exclude 
the possibility that this argument breaks down at the 
crossing point in the exchange path. Related arguments 
can be advanced in which the exchange paths involve 
particle creation at different times, but do not appear to 
us to be conclusive. Our inability to settle this point on 
the basis of simple general arguments has motivated the 
numerical calculations we now explain. 

The identification of the set of states as edge exci- 
tation states of a particular incompressible state in a 
finite-size many-fermion spectrum is both a challenge 
and an important source of uncertainty for the conclu- 
sions we reach. For v = l/{2p-\- 1), the identification is 
accomplishedi^ by appealing to LaughlinSi^ to conclude 
that the low-energy edge-excitation states appear at an- 
gular momenta above Lq = (2p+l)A^(A^— 1)/2, and that 
the dimension of subspaces at fixed angular momentum is 
related to their excess momentum by counting the num- 
ber of modes in a chiral boson Hilbert space^S i.e., one 
state with angular momentum f , two with angular mo- 
mentum 2, three with angular momentum 3, five with an- 
gular momentum 4, et cetera. Previous numerical worki^ 
has verified xLL predictions for i/ = l/{2p-\-l) edges, but 
not at filling factors within the range 1/3 < < 1, where 
experiment and theory appear to be at odds. 

Incompressibilities at many other values of the fill- 
ing factor, e.g., for — m/{2mp± 1) with positive in- 
teger m, have been explained using various hierarchi- 
cal schemea^^i^ and using composite-fermion theoryi2^ 
Composite-fermion theory and variational wavefunc- 
tions based on hierarchy-theory ideas make identical 
predictions'^^ for the values of total angular momentum 
Lq at which maximum- density incompressible states 
(those with edge subsystems in their ground states) ap- 
pear and for the number and chirality of the boson 
branches in their edge-excitation spectra. For practical 
reasons explained more fully below, we limit our study to 
QH systems with two branches of edge excitations that 
have the same chirality, i.e., that propagate along the 
edge in the same direction and have, for our circular 
droplets, excess angular momenta of the same sign. This 
is the case for QH systems at filling factors v = 2j (4p-|-l) 
where the edge spectrum is expected to be that of two 
boson modes that share the same chirality. 

The angular-momentum values at which maximum- 
density states occur depend on the number of electrons 
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TABLE I: Total angular momenta Lo{N, Nqp) for the ground 
state of compact fractional-QH systems as predicted from 
Eq. for p = 1. These states can be regarded as A''- 
electron states in the lowest Landau level that consist of a 
QH droplet at filling factor 1/3 supporting a compact daugh- 
ter droplet of A'^qp quasiparticles. In an equivalent descrip- 
tion, these are A'-composite-fermion states with lowest and 
first-excited composite-fermion Landau levels having occupa- 
tion A' — Aqp and Aqp, respectively. Numbers in bold type 
indicate that the corresponding ground state in our numeri- 
cal spectra was positively identified as a finite-size i/ = 2/5 
QH state." The ground states whose edge-excitation sectors 
are explicitly used here to test the Fermi statistics ansatz are 
indicated by asterisks. 
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these states given by Cappelli et al^^ based on density profiles. 



N, and the number of particles to be transferred from 
the ground to the upper composite-fermion Landau level 
Nqp. This notation is chosen to suggest the analo- 
gous hierarchy-picture description of the same states, 
which makes identical predictions for the set of angular- 
momentum values at which maximum-density states oc- 
cur. In the composite-fermion language, the number of 
particles in the lower composite-fermion Landau level is 
N — Nqp, and maximum- density finite-size v — 2/(4p+l) 
states appear at Lq = {Nqp — 1) * {Nqp — 2)/2 + {N — 
Nqp){N - Nqp -l)/2 + 2pN{N - 1), where the last term 
comes from the Jastrow factor in Jain's variational wave- 
functions and the first term reflects the reduced angu- 
lar momentum of higher Landau- level states*^ Chang- 
ing Nqp and/or N is analogous to a topological or zero- 
mode excitation in a finite-size one-dimensional elec- 
tron system. We denote the N-electron state with 
iVqp quasiparticles and an edge-density sub-system in its 



ground state by | N, Nqp, 0), and its total momentum by 
Lo{N, Nqp). The following relation can be derived from 
the above composite-fermion expression or from hierar- 
chy theory: 

N{N - 1) 

Lo{N,Nqp) = (2p-fl)^^ L + Nqp{Nqp-l)-NNqp. 

(4) 

For the case p = 1, Table shows Lo{N, Nqp) for these 
finite-size maximum- density states. Finite-size numeri- 
cal spectra for small N exhibit non-degenerate ground 
states at these values of L, and a low-energy excitation 
spectrum at small excess L that corresponds to the ex- 
pected chiral two-branch edge when Nqp is close enough 
to « N/2. (A chiral two-branch edge has two low-energy 
states with excess angular momentum equal to 1, five 
with excess angular momentum 2, ten with excess mo- 
mentum 3, et cetera. As an example. Figure ^ shows 
spectra obtained for N — 8 and 9 particles.) We define 
the edge-state Fock-spacc projection V by retaining for 
each N and A'qp only these states. 

The low-energy Fock space is then the direct sum of 
the low-energy Hilbert spaces for different particle num- 
bers N; each of these is in turn the direct sum of or- 
thogonal subspaces labeled by different values of A"qp, 
the number of composite fermions in the first excited 
Landau level. Our understanding of xLL theory is that 
it attempts to describe the physics of QH systems pro- 
jected onto this subspace. It will prove useful to write 
the projected electron creation operator as the sum of 
separate contributions labeled by the change in the num- 
ber of quasiparticles that accompanies a particle addi- 
tion. We therefore write V'^ = Xln ''/'n where the operator 
ipj^ changes Nqp by n. In other words, an electron can 
be added to an iV-particle, TVqp-quasiparticle finite-size 
state in many different ways, distinguished both by the 
resulting edge disturbance and by the number of quasi- 
particles in the resulting (A^-l-l)-particle state. If the cor- 
responding projected electron creation operators tp}^ sat- 
isfy Fermi statistics, their anticommutator {V'l'j V'i} will 
vanish identically. We check for Fermi statistics by eval- 
uating anticommutator matrix elements that involve the 
lowest possible excess momenta and are therefore likely 
to have the smallest finite-size effects. In particular, we 
define 



An,n'{N,Nqp) = {N + 2,Nqp + ji + n',0\c].^rN+i,N,,+nC^K^\N,Nqp,0) , (5a) 
Bn,n'{N,Nqp) = (Ar + 2,iVqp + n + n',0|4^7'w+i,JV,p+«'4j^,^qp,0) , (5b) 

with K, = X^") {N, Nqp) and K2 = k'^"^ (iV -f 1, A^qp + n) where if^"^ {N, Nqp) := Lo{N + 1, Nqp + n)-Lo {N, Nqp) . 
The projector Vn.n^p projects onto the subspace of the iV-particle low-energy subspace having A'qp quasiparticles 
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FIG. 1: Low-energy portion of the exact finite-size spectrum for N = 8 and N = 9 particles with short-range interactions. The 
non-degenerate states at L = 70 (A'^ = 8), L = 87 {N = 9), and L = 92 {N = 9) can be identified as maximum-density ground 
states of a, V — 2/5 finite-size QH droplet. Their energy is separated by a large gap from a continuum of bulk-excited states. 
States at higher angular momenta with energies below that gap (indicated by asterisks) can be unambiguously identified as 
edge excitations. Note that the multiplicity of these states at L = 7f and L = 72 in the left panel and L = 88, 93, and 94 
in the right panel is exactly as expected for a two-branch chiral boson system. For L > 73 in the left panel and L > 89 and 
L > 95 in the right panel, 'edge' states exist that have energies as high as bulk excitations. In the calculations described 
below, we try to select anticommutator matrix elements that do not require us to utilize Fock-space sectors with missing edge 
states. Cases for which there are missing states are explicitly identified in our discussion. The energies plotted here do not 
include the contribution, proportional to total angular momentum, from our model's parabolic confinement potential which 
lifts the energies of states with larger L. The ground state is determined by the strength of this confinement potential. In the 
thermodynamic limit, it will be one of the i/ = 2/5 maximum-density-droplet states when the chemical potential lies in the 
v = 2/5 gap. 



(i.e., composite fermions in their first excited Landau level). It is then obvious that 

An,n' {N, TVqp) = (iV + 2, iVqp + n -f n', | 4j iV + 1, iVqp + n, 0) (iV -f 1, iVqp + n, | 4^ | N, N^p, 0) , 
Bn^n' {N, TVqp) = ^ (iV + 2, iVqp + n + n', I 4j iV + 1, iVqp + n', {t?}) {N + 1, TVqp -t- n', {r,} \ c\^^ \ N, iVqp, 0) 



r 



where | iV" -|- 1, A^qp -I- n', {rf}) are edge excitations in the 
(A^-l-l)-particle, (iVqp-l-n')-quasiparticle system for excess 
total momentum 



(n,n') 



= Ky{N + 1, TVqp + n, n') - KpiN, A^qp, n') , 
= 2p+l + 2nn' - ?i - n' . (7) 



Fermi statistics is satisfied if the ratio B/A equals —1. 

Our numerical results are summarized in Tabled We 
start by considering a QH droplet at filling factor 1/3, 
i.e., states with Nqp = 0. (See the first row in Table ITTI'l 
Although only one case is presented there, it turns out 
that for any finite-size system with short-range inter- 
actions, the electron operator projected onto the sub- 
space of low-energy excitations above the ground state 
at = 1/3 satisfies Fermi statistics exactly. The Fermi 
statistics of edge excitations is a property of Laughlin 
many-body wavefunctions,'^ which are exact in the case 



of the hard-core- interaction model at v = l/(2p -|- l)i^ 
It is apparent from the remaining entries in Table |n] 
that this property of hard-core-model single-branch QH 
edges does not, in general, hold in the f = 2/5 two- 
branch case. We have tested edge-excitation subspaces 
for systems with N = 6, 7, 8, and 9 electrons and find 
that the ratio B/A that 'measures' Fermi statistics often 
differs substantially from —1. Instead, a strong depen- 
dence of B/A on the change in Ngp that accompanies 
particle addition is apparent from our data. For exam- 
ple, consider the sequences i(-^)-L(^+i)-L'^^+^^ where 
particles are added without adding hierarchy quasipar- 
ticles. In composite-fermion language, this corresponds 
to adding electrons to the lowest composite-fermion Lan- 
dau level only. The finite sequences in this class that we 
have studied, e.g., 35-51-70, 51-70-92, and 66-87-111, 
are gathered in the first block of Table ITll Although vio- 
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TABLE II: Results for the ratio between the ma- 
trix elements A = (4^+^' jc|^^P c|^JL<,'^'> and B = 

(Lq'^"''^'' Ic^f^P c|f^ ILq'^'). This ratio should always equal 
— 1 if the edge-electron operator satisfies Fermi statis- 
tics. We have chosen values Ki and if 2 such that 

A = (ir^Vkl4"^^'>(4"+''l4j4"') involves ma- 
trix elements between maximum- density states only. When 
the order of edge-electron operators is reversed, B — 

E.(ir+''l4jit"-'^'>(£5"+''!4j4"'>, with the^sum 
ranging over matrix elements between edge states IL^.^^^') 
with excess angular momentum AK = Lg^' -I- K2 — L^'^'^^K 
The first block corresponds to -Bo,o/^o,o, the second to 
Boa/Ao^i and Bi, 0/^1,0, and the third to Bi^/Aij. 
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"The edge-electron operator for any finite-size Laughlin state sat- 
isfies Fermi statistics perfectly, 
'incomplete set of edge states. 
"Nearly complete set of edge states. 



lations of Fermi statistics are numerically unambiguous, 
they are large only for those cases where the edge sec- 
tor is seriously incomplete (only 5 out of 10 edge states 
could be identified in the L=38 sector for the sequence 
22-35-51). To some extent, this result is not surprising 
since the addition of an electron to the lowest composite- 
fermion Landau level is expected to involve excitations 
of this level only, as in the 1/ = 1/3 case. It seems plau- 
sible that the results should be similar. The deviation 
from Fermi statistics is suspicious, however, when con- 
trasted with the perfect accuracy seen for the = 1/3 
QH droplet. When an anticommutator involving "^j (in- 
creasing the number of higher-Landau-level composite 
fermions by one) is considered, as in the sequences 39- 
51-66, 56-70-87, and 76-92-111 (see the third block in 
Table . the deviations from Fermi statistics are larger. 
Note that some of the same finite-size edge states are 
used here and in the -ipQ case for which Fermi statistics is 
more closely approximated. It is important to recognize. 



however, that the edge sectors used in all these sequences 
are not complete even for the larger systems considered. 
The incompleteness is related to the fact that there are 
only two electrons in the first excited composite-fermion 
Landau level which implies that the daughter droplets are 
substantially smaller than the parent droplets, exacer- 
bating finite-size difficulties. The strongest evidence that 
Fermi-statistics relations are not satisfied comes from the 
remaining sequences in Table ^1 for which the numeri- 
cal edge-state sector is usually complete. There is no 
unique way of estimating a thermodynamic limit for B /A 
from our data. A consistent scheme would have to keep 
N — 2iVqp ^ as iV ^ 00, in order to maintain parent 
and daughter fluids that are similar in size. The small 
system sizes tractable using current numerical methods 
render such an extrapolation impossible. Examining the 
trends in Table however, we are reasonably confident 
that differences |1 -I- B/A\ > 0.05 are significant and not 
merely due to finite-size effects. We note that in the 
case of conventional one- dimensional electron systems, 
a corresponding calculation would always result in exact 
conformation with Fermi statistics; there are no finite- 
size corrections. The v = 1/3 hard-core-model case also 
produces results in agreement with Fermi statistics with- 
out finite-size corrections. In contrast, our calculations 
exhibit large deviations from Fermi statistics for many se- 
quences where the edge states expected for a two-branch 
boson system are clearly resolved in the spectrum, not 
only for those with an incomplete edge sector. This is 
the case for 51-66-87 for example. Furthermore, when 
edge states exist that have energies above the gap for bulk 
excitations and can no longer be clearly identified (this 
is the case, e.g., for 8-particle states at total momentum 
73), the value of B/A turns out to be determined almost 
entirely by contributions from edge states with energies 
below the gap. Inclusion of any number of states (bulk or 
edge) above the gap energy changes B/A only by a few 
percent. Finally, deviations from Fermi statistics do not 
seem to diminish with increasing particle number. On 
the contrary, the more unambiguously identified edge- 
state sectors at larger A^ yield values of B/A that dif- 
fer consistently from —1, especially for anticommutators 
involving the component of the creation operator that 
increases the number of daughter quasiparticles, tpl- 



III. PROPERTIES OF CANDIDATE 
BOSONIZED ELECTRON OPERATORS 

Our numerical results clearly support a multi-branch 
chiral-boson form for the excitation spectrum of a 
fractional-QH edge, but raise new questions about the 
representation of projected electron creation operators 
in terms of these boson fields. This issue is discussed 
in the following section. We start by carefully examin- 
ing the arguments that have been made in x^L theory 
to obtain bosonization identities. Some of these heuris- 
tic arguments must be ruled out if the edge projection 
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of the electron operator does indeed not satisfy Fermi 
statistics. Alternative proposals are discussed, but we 
have not been able to find a simple form that is consis- 
tent with our numerical results, suggesting that the true 
expression may not be universal. 

In a conventional one-dimensional system, the low- 
energy projection of the electron operator -010 
pressed as the sum tpl-^ — V'r + V^l '^^ right-moving and 
left- moving chiral fermion contributions V'rl- ^'^^ these 
chiral fermion operators, an identity relating them to 
the bosonic charge fluctuations of an interacting system 
can be derived rigorouslyi^SiSLSS Fractional-QH edges 
do appear to be realizations of chiral one-dimensional 
systems, as indicated by the multiplicity of low-lying 
many-electron states in our numerically obtained spec- 
tra. Given this observation, one is tempted to push the 
analogy further and search for a bosonization identity for 
the projected edge-electron operators tpt- Assuming that 
it is local in the angular coordinate 9 along the edge, it 
should read 

where z denotes a normalization constant. The 
'Klein factor' is a ladder operator that con- 

nects many-particle ground states, Ul\N,Nqp,0) — 
I -(- 1, iVqp -f n, 0) and commutes with bosonic edge- 
density operators. The correct commutation relations 
for operators with different n, whatever they are, have 
to be encoded in these factors. The chiral phase field 
(j)n{9) is a superposition of edge-density fluctuations. For 
V = 2/{Ap+ 1), the following decomposition of the phase 
field in terms of eigenmodes is always possible: 

0„(0) = ^(A(^)(^?)+en0(")(^) , (9) 

were 0^'^' is the phase field of the charged edge-magneto- 
plasmon mode which corresponds to fluctuations in the 
total edge-charge density, and 0^'^^ is its orthogonal com- 
plement, the so-called neutral mode. The prefactor l/v^ 
of the charged mode in Eq. 10 is mandated by the fact 
that the addition of an electron necessarily increases to- 
tal electric charge by unity. Additional assumptions are 
necessary, however, to fix the values of 

Within xLL theoryjS, the operators ipo and ip2p+i are 
believed to be special in that they create electrons lo- 
calized at the putative 'outer' and 'inner' edges which 
are the boundaries of the outer parent and inner daugh- 
ter QH droplets with filling factors = l/{2p+ 1) and 
— l/[{2p + l)(4p+ 1)], respectively, that comprise the 
= 2/(4p+ 1) QH state. The density fluctuations at the 
'outer' and 'inner' edges are given by {uo/ L) decj)^ and 
(y\l V) dg4>2p+i, respectively, and the definition of charged 
and neutral modes implies that 

(j)^""^ = ^ (i/i02p+i + i^o0o) , (10a) 



^ . (lOb) 



The addition of electrons to the edge with concomitant 
change of 2p -I- 1 — n flux quanta is viewed as adding 
the electron to the 'outer' edge and transferring at the 
same location n fractionally charged quasiparticles from 
the outer QH droplet to the inner one. This suggests 
the relation ijjj^ ~ ■(/'q exp{ini'o{(f>o — <p2p+i)} which is 
equivalent to 

C„ = -^(2n-l) . (11) 

The chain of arguments leading to Eq. (|ll|l involves 
several assumptions that are not obviously satisfied. For 
example, it is not clear why changes in Nqp that accom- 
pany electron addition have to occur by transferring lo- 
calized quasiparticles from an 'inner' edge to an 'outer' 
one. Relaxing this condition would lift the restriction ex- 
pressed by Eq. Hll|l . In fact, two of the present authors 
suggested the different choice = for strongly corre- 
lated fractional-QH edgesi^ The fact that Eq. l(TT)l en- 
sures Fermi statistics has been regarded^ as strong sup- 
port for this line of argument. However, the numerical re- 
sults discussed in Scc.llllsuggest that there is no basis for 
this requirement. Assuming that a simple bosonization 
identity of the form given in Eq. JHJ holds, we can extract 
the values « 0.51 ± 0.02 and £,f « 0.65 ± 0.06 from the 
data for matrix elements of anticommutators {'4'q, V'J} 
and {■(/'I, (See the Appendix for details of the cal- 

culation.) The value of S^f deviates significantly from that 
expected within xLL theory (1/2) because 0| does not 
satisfy Fermi statistics. While the standard deviations 
from the average extracted ^ are reasonably small, we 
have to caution the reader by noting that the assump- 
tion of a simple bosonization identity would imply sym- 
metry of Bn^n' /^n,n' Under exchange n ^ n' . Clearly, no 
such symmetry is exhibited in our data. The significantly 
large deviation between Bq i/Aq i and i?i^o/^i,o raises a 
big question mark: unless rather complex features are as- 
cribed to the Klein factors, a consistent interpretation of 
our data using a local bosonization formula like Eq. (O 
is impossible. 

IV. DISCUSSION AND CONCLUSIONS 

We have shown, by explicit numerical calculation of an- 
ticommutator matrix elements, that the projection of the 
lowest-Landau-level electron creation operator onto the 
low-energy edge-excitation Fock subspace of a = 2/5 in- 
compressible quantum Hall state does not satisfy Fermi 
statistics. We observe a consistent dependence of the 
anticommutation rules on the particular procedure for 
adding electrons, i.e., on the change in quasiparticlc num- 
ber that accompanies particle addition. We find that 
the numerical data cannot be consistently interpreted 
by assuming any simple generalization of conventional 
bosonization identities. In particular, the expression for 
the electron operator solely in terms of rigid edge de- 
formations (magnetoplasmon modes), which two of us 
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argued for previouslj*^ on heuristic grounds is also not 
supported by our calculations. Since our present study 
was performed for a system with short-range interac- 
tions, however, we cannot exclude the possibility that real 
QH samples where long-range Coulomb interactions are 
present may be consistently described by such a bosoniza- 
tion identity, as is suggested by the amazing experimental 
finding that the tunneling-/^ exponent a ~ \/v. 

Our numerical study demonstrates that the specific 
form of the boson representation of the edge electron 
creation operator cannot be inferred by postulating 
Fermi statistics of the projected edge-electron operator. 
Alternatives^^ to the xLL expressions^ cannot be ruled 
out on these general grounds. However, our present data 
for the short-range-interacting case supports neither the 
conventional chiral-Luttinger-liquid picture nor any sim- 
ple alternative. This points strongly towards the pos- 
sibility that there is no simple universal local bosoniza- 
tion identity for the edge-electron operator, a conclusion 
also reached in an independent recent study4i If true, 
this likely implies that electronic properties of fractional- 
quantum-Hall edges depend crucially on sample specifics. 
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APPENDIX: ANTICOMMUTATORS 
CALCULATED FROM BOSONIZATION 



Acknowledgments 

Within the bosonization formalism, the quantities 
Useful and stimulating discussions with A. Auerbach, v4„^„/ and defined by Eqs. (jSJ, do not depend on 

J. T. Chalker, H. Fertig, F. D. M. Haldane, B. I. Halperin, N or TVqp. We find A„_„' = z and 



^ = Xn,n' E <^ + 1' + I I iV -f 1, iVqp -f n', {77}) {N + 1, TVqp + n', {77} | e^^^^'^o) \N + l,N^^ + n' , 0) , 

{r,} 

(A.l) 

where | iV -I- 1, A^qp -t- n', {77}) are excited states with excess (boson) momentum equal to AATp"'"-*. Contributions 
from Klein factors are contained in Xn,n' which satisfies Xn,n = 1- To better understand the two- branch case, it is 
useful to first consider the situation when only a single branch of edge excitations is present. 



1. Single- branch case 



At the edge of a QH sample at filling factor equal to l/(2p -I- 1), a single branch of edge excitations exists which 
corresponds to the edge-magnetoplasmon (charged) mode. [Within the formalism described above, we would have to 
set ^„ = Q,4'n= (j)^^^ = (j>, and let v l/(2p -I- 1) in order to describe the single-branch edge.] Denoting by Qq (oq) 
the bosonic creation (annihilation) operator for an edge-magnetoplasmon excitation with wave number Q, we can give 
expressions for both the phase field entering the bosonization identity and the excited states appearing in Eq. (|A.1|I : 



E 



iQe 



Q>0 



|A^+1,0,(0) = 



JlIiJ^I^ + 1,0,0) 



(A.2a) 
(A.2b) 
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with ''"^r = Ai^Tp = 2p + 1. Here, r > and lr>Q are integers, and (Z) is a partition of the integer Ai^p labeling 
possible many-body excited states. Straightforward calculation yields 

(A^ + 1, 0, I e-''^(0) 1^ + 1,0, (/)) = (A^ + 1, 0, (/) | e"*^'("' | + 1, 0, 0) = [{-vr)^^ . (A.3) 



(/) r-G(0 

^^^Ef^) , (A.4b) 



Inserting this result into Eq. (jA.ljl . we find 

B 



AXp! 

where {l)s denotes a partition with s cycles, and 5^]^^ denotes Stirling numbers of the first kind.^^ (Klein factors 
are irrelevant for calculating B/A of a single-branch edge.) Using relation 24.1.3 from Ref. and introducing the 
generalized binomial coefficients, 

( a \ _ T{a + l) 
\ n 



r(a -t- 1 — n)n! ' 
we find 

Obviously, because v'^ = AKp ~ 2p + 1, we find that Fermi statistic holds for the bosonized expression of the 
projected edge-electron operator at Laughlin-series filling factors: B/A = —1. 



2. Two- branch case 



Lets now consider again the case v — 2/(4p + 1) where two edge branches with the same chirality are present. 
Many-body excited states entering Eq. ljA.l|) are then the eigenstates of the quadratic edge-density-wave Hamiltonian 
which are spanned by bosonic operators that are, in general, some orthogonal linear combinations of the charged and 
neutral-mode excitations. The coefficients in the linear transformation that relate the charged and neutral modes to 
the bosonic normal modes of the Hamiltonian depend on microscopic details, i.e., the velocities of the charged and 
neutral modes as well as their coupling via interactions. However, it turns out that we do not need these coefficients 
explicitly. It suffices to know that the phase fields entering the bosonization formula can be expressed in terms of the 
normal modes, 0„ = '^'j=i '^j"^ where 

ryj"-* — — r= ^„ sin a , (A. 6a) 



r„) _ sma ^ ^ ^ (A. 6b) 

The fields ipj are defined in terms of the bosonic operators aj^Q that annihilate excitations of the jth normal mode 
in analogy to Eq. l|A.2a(l . and microscopic details of the edge determine the value of a. Excited states entering 
Eq. HA.1|I are direct products of excited states in the two normal-mode subspaces whose combined excess momenta 

equals Aifp . Within each of the normal-mode sectors, total excess momentum is partitioned among the possible 



many-boson states as in the single-branch case. We can, therefore, employ the result (|A.5p and find 

R , / (n) («') \ / (") ("') \ 




X«,n' ( '^^ '^^ j . (A.7b) 



Here we have used the addition theorem for general- Ref. 14^ ) Straightforward inspection shows that the sum 
ized binomial coefficients. (See, e.g.. Section 12.2 of 
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Si ^ — ^ +Cn,£,n' is universal, i.e., independent 

of a. Using that, we can write 

^-'^ [ a4-"') ) ■ ^^-'^ 

With the choice of values for ^„ according to xLL the- 
ory, which is given in Eq. one finds + ^„^„' = 
Aifp"'"-*. Since Aifp"'"-* is always odd, and with a 
suitable definition of Klein factors^ such that Xn,n' = 



one obtains B/A = -1. 

Two of the authors suggested previouslySI the choice 
= 0. For the special case of p = 1, our result IjA.Sp 
would predict Sq^oMo.o = —5/16. Addi- 

tional assumptions about Klein factors are then necessary 
for a consistent description of i3o,iMo,i and i3i,oMi,o- 

Having obtained data for i?n,n'/v4„_„', from numerical 
calculations, we used Eq. (|A.8|I to extract the values of 
and check whether a consistent description using a 
simple bosonization identity of the form ||HJ| is possible. 
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